Spatially Resolved Magnetization in Bose-Einstein Condensed BaCuSi2 0e: 

Evidence for Imperfect Frustration 
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In order to understand the nature of the two-dimensional Bose-Einstein condensed (BEC) phase 
in BaCuSi2 0e, we performed detailed 63 Cu and 29 Si NMR above the critical magnetic field, 
H c i = 23.4 T. The two different alternating layers present in the system have very different lo- 
cal magnetizations close to H c \\ one is very weak, and its size and field dependence are highly 
sensitive to the nature of inter-layer coupling. Its precise value could only be determined by "on- 
site" 63 Cu NMR, and the data are fully reproduced by a model of interacting hard-core bosons in 
which the perfect frustration associated to tetragonal symmetry is slightly lifted. 

PACS numbers: 75.10.Jm, 03.75.Nt, 67.80.dk, 76.60.-k 



Weakly coupled antiferromagnetic (AF) spin-1/2 
dimers are currently of great interest, since they offer an 
ideal playground to study the physics of hard core bosons 
on a lattice [1] . All these systems have in common a col- 
lective singlet ground state, separated by an energy gap 
from a band of triplet excitations. For magnetic fields 
H larger than the critical field H c \ closing this gap, the 
dimers acquire a triplet, i.e. a hard-core boson density 
controlled by the bare chemical potential, that is the re- 
duced magnetic field h = H — H r \. Depending on the 
geometry of the inter-dimer couplings, a rich variety of 
ground states can occur: canted XY AF order which can 
be described as Bose-Einstein condensation (BEC) [2, 3], 
observed in TlCuCl 3 [3, 4], BaCuSi 2 6 [5], DTN [6], and 
BPCB [7]; magnetization plateaus described as Mott- 
insulators, observed in SrCu2(B03)2 [8]; Bose glasses ob- 
served in doped DTN [9] and other systems [10]; and 
even supersolids [11] for which a physical realization has 
not been found yet. The hallmark of a BEC system is 
the field dependence of the temperature Tbec of tran- 
sition into the ordered phase of canted XY AF type, 
Tbec °c h?, where (f> = 2/d = 2/3 in the standard 3D 
case. Particularly interesting is the case of quasi-2D ar- 
rays of dimers, since formally the BEC in 2D can only 
occur at T = 0. However, depending on the strength 
and the nature of the inter-layer coupling which stabilizes 
the 3D order, the Tbec vs. h dependence can be dom- 
inated by the properties of the Quantum Critical Point 
at T = = h, and exhibit 2D exponent = 1. 

During the last few years BaCuSi2C>6, also known as 
Han purple, has emerged as an excellent candidate for 
studying this quasi-2D case. Its Tbec w &s found to 
vary linearly with h, i.e. <f> = 1 [12, 13], pointing to a 



2D BEC. This was originally explained in terms of per- 
fect frustration between the Bose-Einstein condensates 
of adjacent layers, due to the fact that the Cu dimers 
form a body-centered tetragonal (bet) lattice [14] (in- 
set of Fig. 1). However, it is now well established that 
BaCuSi206 undergoes around 90 K a first order struc- 
tural transition [15], introducing an alternation of dif- 
ferent A and B layers along the c-axis, observed both 
by inelastic neutron scattering [16] and by NMR [13], as 
well as some structural incommensurability. Since per- 
fect frustration relies on the AF coupling of one spin to 
four equivalent neighbors in the adjacent layer, the pres- 
ence of incommensurability is in general incompatible 
with perfect frustration. As a result, theories with differ- 
ent levels of frustration have been put forward [17, 18], 
and, to reach a consistent explanation of the remarkable 
properties of BaCuSi2C>6, the next step is to get reliable 
information on the actual level of frustration. 

In this letter we show that the complexity brought 
about in BaCuSi206 by the presence of two different lay- 
ers provides a unique opportunity to quantitatively de- 
termine and model the extent of interlayer frustration. 
Indeed, as shown in Refs. 17, 18, the magnetization, i.e. 
the boson density, of the less polarized layers just above 
H c i depends crucially on the level of frustration: it is ex- 
pected to grow quadratically with h if the frustration is 
perfect [18], and linearly otherwise [17]. In addition, the 
bosons are expected to remain uncondensed close to the 
first critical field if frustration is perfect [19], while con- 
densation will occur immediately otherwise [17]. So, the 
only possibility to directly distinguish between the differ- 
ent models is a precise determination of the local mag- 
netization of the less polarized layers. In that respect, 
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the 29 Si NMR data of Rcf. 13, which have unambiguously 
demonstrated the existence of two different contributions 
to the spectrum, hence of two types of layers, are not suf- 
ficiently precise for our present purpose, mostly because 
the lines of the two types of layers overlap. 

To overcome this difficulty, we have decided to per- 
form low temperature 63 Cu NMR measurements in the 
vicinity of the quantum critical field H c i = 23.4 T of 
BaCuSi 2 6 [13]. Unlike the 29 Si NMR spectra [13], the 
63 Cu NMR signal from the B layers is clearly separated 
from the one from the A layers, and provides an accu- 
rate determination of boson densities in the B layers as a 
function of h (Fig. 1). As predicted in [18], a quadratic h 
dependence has been observed, however, its size is much 
stronger than the predicted one, and there is also a non- 
negligible linear component. We attribute this latter to 
a deviation from the perfect frustration due to the in- 
commensurate (IC) lattice modulation [15], and a weak 
breaking of the tetragonal symmetry [20] , both of which 
characterize the structural distortion at 90 K. 

BaCuSi206 can be described as Cu 2+ spin- 1/2 dimers 
located on a bet lattice with Cu-Cu bonds perpendic- 
ular to the a6-plane (inset of Fig. 1). The real space 
group at room temperature is indeed IA\/acd [21]. How- 
ever, below 90 K a first order structural transition ren- 
ders adjacent layers inequivalent [20], and also induces 
an IC distortion of the lattice along the 6-axis [16]. This 
IC distortion converts the typical splitting of each NMR 
line below Tbec [7] into a continuously distributed spec- 
trum [13]. In addition, although the boson densities ua 
and nB in the A and B layers are very different, the 29 Si 
NMR lines originating from the two types of layers over- 
lap due to the small sensitivity of the 29 Si NMR to local 
boson density. Finally, through a significant long-range 
(dipolar) coupling, the 29 Si nuclei from the B layers are 
also somewhat sensitive to the A layer boson density. Al- 
together, this renders an accurate determination of the 
small n-Q values by 29 Si NMR quite difficult. Although 
experimentally more complicated, a study by 63 Cu NMR 
presents major advantages. The 63 Cu hyperfine tensor is 
much larger, 63 A ZZ = —16.4 T//xb [13] (defined later, see 
Eq. 2), enhancing the sensitivity of 63 Cu NMR line posi- 
tion to small values and variations of local boson density. 
One thus expects that above H c \ the different «a and 
7iB values correspond to separated NMR lines, rendering 
the determination of ub much easier. 

The 63 Cu and 29 Si NMR have been performed in the 
23-27 T field range, in a resistive magnet at LNCMI. 
The sample was a 2. 4x2. 4x1. 2mm 3 single crystal of 
BaCuSi2 0g, grown in a mirror furnace, and fully enriched 
in 29 Si to enhance the 29 Si NMR signal. It was placed 
in the mixing chamber of a 3 He- 4 He dilution refrigerator, 
with its c-axis oriented parallel to the external field H, 
i.e., c || z. NMR spectra have been recorded at fixed field 
by summing up Fourier transforms of the NMR echoes 
taken at regular frequency steps. The frequency v for 



the 63 Cu NMR spectra corresponds to the central line of 
63 Cu, and is taken relative to the reference uq = 63 j e ffH, 
where 63 ~/ cS = 63 ~f(l + K olh ), with 63 7 = 11.285 MHz/T 
and iforb = 1-68 % [13]. The value of the applied field H 
was calibrated with the metallic 27 A1 reference placed in 
the same coil as the sample. 

In a Bosc-Einstein condensate spins (S) bear both a 
longitudinal M z = g z ^B(S z ) and a transverse, staggered 
component of the magnetization M± = g±fJ-B{S±) , where 
g z and g± are the corresponding components of the g 
tensor. In general, both components define the observed 
NMR frequency through the corresponding components 
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FIG. 1: (color online) Magnetic field dependence of the low- 
temperature (0.05-0.08 K) 63 Cu NMR spectra on entering the 
BEC phase of BaCuSi2 06 (green lines). The vertical offset of 
the spectra defines the corresponding field value. The spectra 
are almost exclusively due to B layers, and their first moment 
(open circles) measures the corresponding local magnetiza- 
tion. The dotted line is the theoretical prediction, see the text 
and Fig. 2. The x5 zoom of selected spectra (orange lines) 
reveals a weak, broad and strongly shifted contribution of A 
layers, and the orange triangles denote a rough estimate of its 
first moment. Inset shows schematic structure of the spin 1/2 
dimers: green lines present perfectly frustrated, symmetric 
inter-layer exchange couplings, and dotted lines denote weak 
breaking of this symmetry. 
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of the hyperfine coupling tensor A zz and A z ±: 

v± = vv+ jA zz M z ± jA zA _M ± . (1) 

While for the 29 Si NMR in BaCuSi 2 6 this is indeed 
the case, the hyperfine tensor of the Cu is dominated by 
on-site interactions (core polarization and on-site dipolar 
interaction) and is diagonal with its main component A zz 
along the c-axis. The Cu line shapes are thus insensitive 
to the staggered magnetization (as long as H \\ c) and the 
frequency shift of the first moment v\ (average frequency) 
of the line, Av = Vi — vq, reflects directly the longitudinal 
magnetization only, that is the boson density ^a(B) by: 

n A(B) = 2(S^} = 2A t , A(B )/( 63 7.g cc MB 63 A ZZ ), (2) 

where g cc = 2.30 [23]. This formula implies that Uacb) = 
1 when there is one boson per dimcr in the A(B) layer. 
With this convention, one has ha + tib = 2 at the sat- 
uration field H C 2 « 50 T. Fig. f presents the field de- 
pendence of the 63 Cu NMR line in the BEC phase. The 
fact that the boson density as well as its fluctuations are 
much bigger in the A layer makes their contribution to 
the spectra nearly unobservablc for two reasons. First, 
the intensity of the A line is much smaller, due to a much 
shorter spin-spin relaxation time T 2 reflecting the fluctu- 
ations. This mechanism is particularly effective close to 
the phase transition. Second, the width of the line, being 
related to boson density, is expected to be much broader 
for the A line and thus of proportionally smaller intensity 
(as the total signal is given by the integral over the line) . 
Indeed, in the enhanced scale spectra in Fig. 1, the A 
line can be observed as a very weak broad signal on the 
left of the main B line only above 23.9 T where the two 
lines become fully separated. Note that the main line 
corresponding to the B layer is obviously shifting much 
more slowly than the one corresponding to the A layer, 
reflecting a much lower boson density Ub- 

Using Eq. (2) , the first moments of the lines shown in 
Fig. 1 are converted to «b and their field dependence is 
plotted in Fig. 2. As ha determined from 63 Cu NMR 
suffers from huge error bars due to a very poor signal to 
noise ratio of the A lines, we have instead plotted the 
total boson density n A + «b , deduced from the first mo- 
ment of the 29 Si line-shapes recorded on the same sam- 
ple [22]. As already mentioned, the contributions of the 
A and B layers overlap in the 29 Si spectrum, and the 
lincshape is sensitive to both the transverse magneti- 
zation and IC distortion [13]. Fortunately, both latter 
effects average to zero when calculating the first mo- 
ment V\. Therefore, the 29 ^i of the total Si spectrum, 
including the A and B layer contributions, gives an ac- 
curate determination of n A + through the formula: 
"A + n-B = 4( 29 z^i - 29 i/ )/( 29 75ccMb 29 A ZZ ), which is the 
same formula as Eq. (2) corrected for a factor of 2 for 
proper normalization. A possible source of error is the 
difference in the 29 A ZZ values for the A and B layer; here 
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FIG. 2: (color online) Magnetic field dependence of boson 
densities. The tlb from 63 Cu NMR (open blue circles) and 
tia + nB from 29 Si NMR (red squares) are compared to three 
different theoretical models (see the text) having perfect, par- 
tial and no frustration (thin black, thick green, and dashed 
magenta lines, respectively). Thin cyan dashed line shows 
the boson density of identical, decoupled layers (2D). Inset: 
Zoom on hb in the vicinity of H c \. 



we have used the former one (0.30 T / Ub) to ensure that 
the error is negligible in the relevant field range where 
n B < n A [22]. 

We now focus on the variation of tib as a function of 
h. In earlier models [14] all the layers were taken to be 
identical, which is not realized in BaCuSi20g. When the 
finite energy barrier, that is the difference between the 
two local gap values, A = Ab — Aa ~ 0.4 meV [13, 16] is 
considered, perfect frustration along the c-axis leads to a 
classical decoupling of the layers, and to the absence of 
bosons in B layers close to H c \ [17], a result in contradic- 
tion with our previous [13] and new 29 Si NMR data [22], 
the latter ones being more accurate since recorded on a 
29 Si enriched sample. However, as shown in Refs. 18, 19, 
quantum fluctuations (spin-wave corrections) above such 
a classical ground-state induce above H c \ a small boson 
density tib oc h 2 . As a corollary, a full three dimensional 
coherence is restored in this case through an effective 
inter-layer hopping of bosons, cx h, leading to a mod- 
ified exponent Tbec tx h^ =1 . in agreement with experi- 
ments [12, 13]. 

In the inset of Fig. 2 we plot in detail tib (ft-), which 
is an order of magnitude smaller than ua , and whose /in- 
dependence obviously exhibits a quadratic contribution. 
However, the measured density ub turns out to be much 
larger than earlier theoretical predictions based on a per- 
fectly frustrated scenario [18, 19]. Moreover, a linear 
component is also present in the observed «b(/i) depen- 
dence. As already mentioned, the key to explain these 
deviations is an imperfect frustration between A and B 
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layers, which has two origins. According to the recent 
low-temperature average structure of BaCuSi20g [20], 
the tetragonal symmetry is broken, which introduces two 
slightly different Cua-Cub distances, corresponding to 
the two type of Cua sites having opposite staggered mag- 
netization. This provides a direct and site-independent 
lifting of the perfect frustration. Another source of im- 
perfect frustration is the IC modulation in the b direction, 
which is along the diagonal connecting the next nearest 
neighboring dimers [15]. While both effects are very dif- 
ficult to estimate, the NMR data provide, with the help 
of the following theoretical analysis, a reliable estimate 
of their total "effective" size. 

As discussed previously in Refs. 5, 14, 18, 19, one can 
map the system of coupled spin- 1/2 dimers in BaCuSi20g 
onto a lattice model of hard-core bosons whose quantum 
dynamics is governed in each layer £ = A, B by the Hamil- 
tonian 

U\ = t,| [ b l( r )W) + h -c- + ^(r)^(r') 

n 

- fnJ2n £ (r), (3) 

r 

where the sum runs over the nearest neighbors on a 
square lattice. It is important to notice the positive 
(antifcrromagnetic) hopping term tu and the distinct 
chemical potentials on different layers: [ia = Mb + A. 
The inter- layer (or "transverse") tunnelling term, which 
mixes bosons from layers A and B, reads: 

H± = [fe A (r)fe B (r')+h.c. + nA(r)n B (r') 

X 

+ ^[&i(r)fc B (r') + h.c. + nA (r)«B(r')] ,(4) 

X 

where the first term accounts for the perfectly frustrated 
hopping between A and B, and the second term mimic a 
direct (unfrustrated) tunneling between neighboring lay- 
ers. One can immediately rule out the scenario of identi- 
cal and perfectly decoupled layers (A = t± = t' ± = 0) [14] 
for which the boson density is plotted in Fig. 2 by thin 
cyan dashed line. Although this gives an approximately 
correct description of the total measured density, in par- 
ticular close to H c i , the B layers should simply not exist 
in this case. A realistic model should include different 
layers A and B (A ^ 0) and an inter-layer couplings 
between them. In order to get the best description of 
NMR data, we fix the parameters tn = 0.3 meV and 
A = 0.3 meV, and allow tj_ and t'j_ to vary. A spin-wave 
analysis of this effective hard-core bosonic Hamiltonian 
is performed at the linear 1/5 order [18, 24]. We have 
examined three possible scenarios: (i) perfect frustration 
with t±/tu = 1/8 and t' ± = 0; (ii) no frustration with 
t± = and t' ± /t\\ = 1/20; (hi) partial frustration with 



t±/t\\ = 1/8 and ij_A|| = 1/50. In all these cases, the 
inter-layer barrier A = 0.3 meV has been chosen in or- 
der to give the best description of the densities measured 
by NMR. This value appears somewhat smaller than the 
experimental ones [13, 16], but we expect that the IC 
modulations, particularly strong in B layers, naturally 
lead to a reduction of this internal barrier. We now ex- 
amine in details the three aforementioned scenarios: 

(i) While a perfect frustration is particularly interesting 
since it would induce an exotic boson density in the B 
layers [19], it is unfortunately unable to explain the mag- 
nitude of the density n-Q. This is evidenced in the inset 
of Fig. 2 where we see that the quadratic contribution 
coming from quantum fluctuations is extremely small. 
Moreover, the perfect frustration scenario causes a clear 
cusp in the field dependence of the total density n A + ^b 
around H w 26 T, while in the experimental data, one 
observes only a slight bending. 

(ii) The completely unfrustrated case requires a more de- 
tailed discussion. The present analysis focuses on the 
quantum critical regime where the bosons are very di- 
luted, but we have to keep in mind the full experi- 
mental phase diagram, and in particular the maximum 
?bec ~ 3.8 K value occurring at H « 37 T [5], which 
imposes a transverse tunneling t' ± = 0.05 t\\. This value 
of the direct (unfrustrated) hopping leads to theoretical 
densities tia and tib much larger than the experimental 
data. In particular, tib is about five times bigger, and 
the bending of the n\ + tib curve has completely disap- 
peared with such a scenario, as seen in Fig. 2. 

(iii) We are then left with a third scenario, which incor- 
porates a frustrated coupling t±/tu =1/8 and a much 
smaller direct hopping t' ± /t\\ = 1/50 lifting the perfect 
frustration by 2t' ± / (4t± + 2t' ± ) w 7%. Within such a par- 
tially frustrated model, one can qualitatively and quan- 
titatively explain all the features observed in NMR ex- 
periments: the small but finite density having a pro- 
nounced convexity, and the total density nx + ?ib which 
smoothly departs from a 2D behavior above 25 T. Re- 
garding the scaling of TbecC 1 ), we expect the contri- 
bution of such a small direct coupling to restore a true 
3D critical regime with an exponent <j> = 2/3 only below 
(t'j_) 2 /A fa 10 mK, which is out of the current experimen- 
tal range of measured Tbec(^) in BaCuSi206 [12, 13]. 

To summarize, taking advantage of the presence of two 
types of layers, the combination of 63 Cu and 29 Si NMR 
has been shown to be an extremely sensitive probe of 
frustration in BaCuSi206. In particular, the density of 
the less polarized layers in the BEC phase, revealed by 
63 Cu NMR, is not purely quadratic in (H — H c i) and is 
much larger than predicted assuming perfect frustration. 
This unambiguously demonstrates the presence of imper- 
fect frustration of the inter-layer couplings, and we have 
quantified this effect in a theoretical model, in which the 
unfrustrated component amounts to « 7 % of the total 
inter-layer coupling strength. We hope that the present 
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results will motivate further structural and ab-initio in- 
vestigations to determine the precise form of the mag- 
netic model that describes this remarkably rich and in- 
teresting system. 
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